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Abstract 



The best approximation by bounded product functions is calcu- 
lated for some very simple two- valued functions of two variables. 

Introduction 

Here is a finite-dimensional extremal problem: given an n x n-matrix 



of numbers djj = ±1, maximize Ylij a i,j^i c j over an ^-vectors (foj)j, (cj)j of 
numbers bi — ±1, Cj — ±1. 

The corresponding infinite-dimensional problem is: given a measur- 
able function of two variables / : X x Y — > { — 1,1}, maximize 
JJ f(x,y)g(x)h(y)dxdy over all measurable functions g : X — > { — 1,1}, 
ft, : y — >{— 1,1}. Here X, y are given measure spaces of finite measure. 

Equivalently, we seek the best /^-approximation of a given function /(•,•) 
by factorizable functions g(-)h(-) (all values being ±1). 

More generally, we may consider measurable vector-functions g : X — > M d , 
/i : y — > lR d (for a given dimension d) satisfying \g(x)\ = 1, = 1 for all 

x, y (| • | stands for the Euclidean norm). In this case g(x)h(y) is interpreted 
as the inner product; that is, we maximize 

f(x,y)(g(x),h(y))dxdy. 

(Still, / : X xY — >• { — 1, 1}.) The case <i = 1 is just the scalar case considered 
above. We may also use an infinite-dimensional separable Hilbert space H 
in place of M. d (the case d = oo). 

About a relation to Bell- type inequalities see for example [I], especially 
(2.18), (2.20) and (2.24)-(2.25). The case d = 1 is related to classical Bell- 
type inequalities; d = oo — to quantum Bell-type inequalities in general; 
d = 2 — to quantum Bell-type inequalities for a maximally entangled pair 
of qubits (the singlet state of two spin- 1/2 particles). 
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Recently some physicists [2] got especially interested in two examples (see 
below) of such extremal problems. They conjectured the optimal functions 
g, h (in both examples, for all d) but did not prove optimality of these func- 
tions. For d = 1 they deduce their claims from Inequality (7) of [3] (for 
n — > oo). However, looking at [3j I did not find a proof of (7); it is just 
checked by inspection for some small n. My goal is to prove the optimality. 

Example 1. X — Y — (0,1); 




f{x,y) 



/(-) 
the given function 



1 if x < y, 
-1 otherwise. 



+ 



g(-)h(-)iord=l 
one of the optima 



Example 2. X = Y = T = R/Z (the circle of length 1); 




f{x,y) 



/(-,-) 

the given function 



1 if dist(x,y) < 0.25, 
-1 otherwise. 



g(-)h(-) for d= 1 
one of the optima 



The following result holds for Example [T] as well as for Example |2j 
Theorem 1. 

\ for d — 1, 
2 for 2 < d < oo. 



max // f(x,y)(g(x),h(y))dxdy 



(The maximum is reached for every d.) 



1 Some necessary conditions of optimality 

Let functions g,h maximize If(g,h) = jj f(x,y)(g(x),h(y))dxdy over all 
measurable g : X — > M. d , h : Y — > M d (for a given dimension d) satisfying 
|^(x)| = 1, \h(y)\ = 1 for all x, y. We introduce G : X -> M d , H : F -> M d by 

= y f{x,y)Hv) d V' 
H{y) = / f(x,y)g(x)dx . 
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Clearly, f(g(x), G(x)) dx = I f (g, h) = J{H(y), h{y)) dy. The optimality of g 
(for the given h) implies that g(x) = G(x)/\G(x)\ for almost all x such that 
G(x) ^ 0. Thus, I f (g,h) = J\G(x)\dx. Similarly, h(y) = H(y)/\H(y)\ for 
almost all y such that H(y) ^ 0, and //(<?, h) = f \H(y) \ dy. 

We turn to Example [U X — Y = (0, 1) and f(x,y) = sgn(y — x) a.e. 
(almost everywhere). We have G(x) = — J* h(y) dy + f*h(y)dy, therefore 
G is absolutely continuous and G\x) = —2h(x) a.e. Also, G(0) + G(l) = 0. 
Similarly, H'(y) = 2g(y) a.e., and H(0) + H(l) = 0. We get a system of 
differential equations: 

G'(x) = -2 H(X) 



H'(x) 



\H(x) 
Gjx) 
\G(x)\ ■ 



Each equation holds almost everywhere, except for the points where its right- 
hand side is undefined. 

1.2 Remark. It can be deduced that |G(x)| + |-ff(x)| = const and therefore 
|G(x)| + \H(x)\ = 2If(g, h) for all x G (0, 1). However, this fact will not be 
used. 

Now we turn to Example [2j X = Y = T and f(x,y) = sgn(0.25 — 

dist(x,y)) a.e. We have G(x) = Kv) d V ~ Go.S Kv) d 2/' therefore 

G is absolutely continuous and G'(x) = 2h(x + 0.25) — 2h(x — 0.25) a.e. 
(since x + 0.75 = x — 0.25 in T). Also, G(x + 0.5) = —G(x). Similarly, 
H(x + 0.5) = -H{x). We have G'(x) = AH{x + Q.2h)/\H{x + 0.25) | and 
H'(x) = 4G(x + 0.25)/|G(a; + 0.25)|. Introducing H(x) = H(x + 0.25) we get 
a system of differential equations: 

G'{x) = aM {x) 



;i.3) 

H'(x) = -4 



H(x)\ 
G(x) 



\G(x)\ • 



Each equation holds almost everywhere, except for the points where its right- 
hand side is undefined. 



2 Dimension one 

Let d=l. 
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Differential equations fll.lH describe a dynamics on the plane M 2 ; the 
integral curves evidently are the squares \G\ + \H\ = c, c G [0, oo). 

H H 

/ \ 
/ \ 

► G 

\ / 
\ / 



The solutions are periodic: G(x + 2c) = G(x), H(x + 2c) = H{x). Also, 
G(x + c) = -G{x), H(x + c) = -H(x). The condition G(0) + G(l) = 0, 
H(0) + i?(l) = selects a sequence of solutions: c G {1, |, |, . . . } U {0}. We 
have to choose c as to maximize h). Using the equality \G(x)\ dx = 

I f (g,h) = /„ |/T(ar)| dx we get I f (g,h) = 0.5 ^(^(a;)! + \H(x)\) dx = 0.5c. 
The maximizer is c = 1 and the maximum is If(g, h) = 0.5. 

The starting point (G(0),H(0)) can be chosen arbitrarily on the square 
|G(0)| + |tf (0)| = 1. 

H H 



For instance, choosing G(0) = 1, H(0) = we get > on (0, 1), G(-) > 
on (0,0.5) and G(-) < on (0.5,1). Thus, fc(-) = 1 on (0,1), ^(-) = 1 on 
(0,0.5) and g(-) = —1 on (0.5, 1). This is the solution shown on the picture 
(see Example CP . 

Differential equations ( 11. 3p are quite similar: \G\ + \ H\ = c; the condition 
G(x + 0.5) = -G(x), H(x + 0.5) = -H(x) selects c G {1, |, ~, . . . } U {0} 
again. The maximizer is c = 1 and the maximum is If(g, h) = 0.5. 

H 
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Choosing G(0) = 1, H(0) = we get G(-) > on (-0.25, 0.25), G(-) < on 
(0.25,0.75), H(-) < on (0,0.5) and H(-) > on (-0.5,0). Thus, #(•) = 1 
on (-0.25,0.25), g(-) = -1 on (0.25,0.75), h(-) = 1 on (-0.25,0.25) and 
h(-) = — 1 on (0.25,0.75). This is the solution shown on the picture (see 
Example [2]) . 



3 Higher dimensions 

We start with Example [2] and use Fourier series: G, H : T — > M. d , 
G(x) = a k e 2nikx , H{x) = ^ b k e 2 ™ kx ; 

k k 

a k = f G{x)e- 27T[kx dx , b k = [ H{x)e~ 2 ^ kx - 
Jo Jo 

a k , b k G C d for k E Z. We know that G(x + 0.5) = -G(x) and H(x + 0.5) = 
—H(x), therefore a k = 0, b k = for all even /c. Especially, a = 0, 6 = 0. 
Taking into account that 

/•i 

|2 



/ 1 |G(x)| 2 dx = J] 
/ |G'(x)| 2 dx= (2tt) 2 Vfc 



2 I« I 2 
Oft | j 



we get 

\G ; (x)\ 2 dx>4n 2 I \G(x)\ 2 dx. 



According to (11.31) . |G"(x)| = 4 for almost all x, therefore Jjj 1 \G(x)\ 2 dx < 
A/tt 2 and 



f\G{ 
Jo 

Using the equality If(g, h) = J \G(x) \ dx we get 



2 

x) | dx < — . 

71 



if(g,h)<-. 

This bound can be reached already for d = 2. Namely, we may take two- 
dimensional vector-functions 

(2 2 \ / 2 2 \ 
G(x) = ( — cos 2nx, sin 2nx ) , H(x) = ( sin 2nx, cos 2ttx ) 
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satisfying (11.31) and the conditions G(x + 0.5) = —G(x), H(x + 0.5) = —H(x). 
Then |G(a;)| = 2/tt for all x, therefore If(g, h) = -. More explicitly, 

g(x) = (cos 2ttx, — sin 2ttx) , h(x) = (— sin 2ttx, — cos 2ttx) . 

We turn to Example GQ here G, H : [0, 1] -> G(0) + G(l) = 0, #(0) + 
#(1) = 0. We extend G, if to [0, 2] letting 

G(l + x) = -G(x) , + x) = -H(x) for a; G [0,1] 

and use Fourier series 

g{x) = a ^ ikx , H ( x ) = bke7Tikx 

k k 

where = 0, = for all even k. We have 



2 1 |2 



|G(x)| 2 dx = ^ K| 2 , I |G'(x)| 2 dx = n 2 Y,k 
Jo k Jo k 

4=/ |G'(x)| 2 dx > 7T 2 I |G(x)| 2 dx, / |G(x)|dx<-. 
Jo Jo Jo 71 

The bound is reached (in particular) for the two-dimensional vector-functions 

(2 2 \ (2 2 \ 
G(x) = ( — cos ttx, — sin ttx ) , H(x) = ( — sin ttx, cos ttx ) , 

\7T TT / \TT TV J 

which means 

g(x) = (cos7rx, sin ttx) , h(x) = (sin7ra, — cos7rx) . 
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